THE AGE-DEPENDENCE OF THE DETECTABILITY OF 
COMETS ORBITING SOLAR-TYPE STARS 



M. Jura 

Department of Physics and Astronomy, University of California, Los Angeles CA 

V^j-"' 90095-1562; jura@clotho.astro.ucla.edu 

O 
O 

<N ■ ABSTRACT 
-i— > 
O 

The outflows from comets in orbit around G-type main-sequence stars can be 
detected when they produce transient OH absorption lines in the spectrum near 
3100 A of the host star. There is only about a 3 x 10~ 8 probability of detecting 
an analog to comet Hale-Bopp orbiting an analog to the Sun. However, for young 
solar-type stars with very large numbers of comets, possibly delivering water to 
terrestrial planets, there is as much as a 1% chance that any sufficiently sensitive, 
randomly-timed observation may detect such transient absorption. 
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Subject headings: astrobiology - comets: general - planetary systems 



INTRODUCTION 



A fundamental goal in modern astronomy is to learn more about the origin and evolution 
of the solar system and other planetary environments where life may exist. Circumstellar 
dust around main sequence stars was discovered with IRAS (Aumann et al. 1984) and has 
been studied extensively to learn indirectly about the origin and evolution of larger parent 
bodies which might resemble asteroids or comets (Lagrange, Backman & Artymowicz 2000; 
Zuckerman 2001). To date, however, very little is directly known about these parent bodies. 
The goal of this paper is to describe an observational strategy to identify and characterize 
individual extrasolar comets. 

Water-rich asteroids and comets are probably responsible for delivering to the Earth 
most of its water (Morbidelli et al. 2000). By studying comets and related objects around 
young stars, it may be possible to learn more about this process which is fundamental for 
the development of life as we know it. 



There have been a number of previous discussions about possible observational signa- 
tures of comets around other stars. Alcock, Fristrom & Siegelman (1986) suggested that the 
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presence of metals in the atmospheres of white dwarfs could be explained by the collision of 
comets with the star, but this hypothesis is not yet supported by additional evidence (see 
Zuckerman et al. 2003). Melnick et al. (2001) detected gas-phase H 2 in the outflow from 
IRC+10216, a carbon-rich mass-losing red giant, which they interpreted as being produced 
by the sublimation of comets (Stern, Shull k Brandt 1990, Ford k Neufeld 2001). However, 
the observed water might result from chemical reactions in the dense outflow (Willacy 2004). 
Jura (2004) has argued that the lack of excess 25 fxm radiation in IRAS data for first ascent 
red giants means that these stars typically have less than 0.1 M e of comet-like Kuiper Belt 
Objects in orbits at ~ 45 AU. 

There is compelling evidence for Falling Evaporating Bodies (FEB's) around (3 Pic, a 12 
Myr-old star (Barrado y Navascues et al. 1999, Zuckerman et al. 2001) with a substantial 
amount of circumstellar dust (Artymowicz 1997, Vidal-Madjar, Lecavelier des Etangs k 
Ferlet 1998). Transient absorption lines which can be attributed to infalling planetesimal 
also have been detected around other young A-type stars (see, for example, Grady et al. 
1996, Roberge et al. 2002). Although Lecavelier des Etangs, Vidal-Madjar k Ferlet (1996) 
and Li k Greenberg (1998) have proposed that much of the dust around (3 Pic arises from 
comets, the observed FEB's contain a substantial amount of refractory material, and their 
degree of resemblance to ice-rich comets in the solar system is uncertain (Karmann, Beust 
k Klinger 2001, 2003, Thebault, Augereau k Beust 2003). 

There are very substantial uncertainties regarding the Oort belt comets in the solar 
system. Estimates of the total mass in this region range from 0.6 M e (Stern k Weissman 
2001) to 40 M e Weissman (1996), while there are a factor of 100 fewer observed "dormant" 
comets than predicted with the simplest models (Levison et al. 2002). 

Given the major unknowns, in this paper we describe schematic rather than detailed 
models. We consider the production of gas-phase absorptions in the observed spectra of main 
sequence stars caused by transiting comets which are warm enough that water-ice sublimates 
(see Beust, Karmann k Lagrange 2001). Previously, the effects of the photometric variations 
by cometary dust have been computed by Lamers, Lecavelier des Etangs k Ferlet (1997) 
and Lecavelier des Etangs, Vidal-Madjar k Ferlet (1999); here, we consider the possibility 
of detecting comets by their spectral line absorptions (see also Smith, Black k Oppenheimer 
1981). We show results for OH because it is abundantly produced in comets and because it 
can be detected from the ground. 
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2. MODEL 

2.1. Outflow from the Comet 

Following Whipple & Huebner (1976) and Festou (1981), we adopt a simple model 
for the gaseous outflow from comets. When a comet approaches its host star, the ice is 
sublimated, and subsequently the gas-phase H 2 is photodissociated mainly into OH and H: 

hu + H 2 + OH (1) 

The resulting OH is mainly photodissociated into O and H: 

hv + OH^H + O (2) 

Finally, the O and H atoms are ionized, typically by charge exchange with the wind particles 
from the host star. Below, we compute the spatial distribution of the OH molecules. 

We assume that the comet has an active area, a, and lies at a distance, R, from the 
host star. [We use capital and small letters to denote locations and speeds in the coordinate 
system centered on the host star and comet, respectively] We assume that the comet has 
a single surface temperature, and we focus on the parameter regime where the comet is 
sufficiently hot that cooling by ice sublimation balances heating by radiation from the host 
star and radiative cooling from the comet is relatively unimportant. If AE is the average 
energy required to remove a water molecule from the comet, and denotes the luminosity 
of the host star, then the rate of production of H 2 molecules, Nh 2 o, is given by: 

The sublimation energy of pure water ice at K is 7.6 x 10~ 13 erg/molecule (for example, 
Ford & Neufeld 2001). However, not all of the incident stellar energy is converted into ice 
removal. We adopt AE = 2.0 x 10~ 12 erg/molecule, approximately the result given by 
Sekanina (2002) for the erosion energy of Sun-grazing comets. With equation (3) and our 
adopted value of AE, we reproduce within a factor of 2, the rates of H 2 and OH production 
given by A'Hearn et al. (1995) and Crovisier et al. (2002) for ~100 comets studied since 
1976. For reference, the comet with the largest measured outflow was Hale-Bopp where the 
water molecule production rate at 1 AU was approximately 1.0 x 10 31 s _1 (Dello Russo et 
al. 2000, Makinen et al. 2001, Harris et al. 2002, Morgenthaler et al. 2001 and references), 
thus implying a = 1.5 x 10 13 cm 2 . 

To use equation (3), we assume that the comet is warm enough that sublimation is 
the dominant cooling process. This assumption fails for comet temperatures below ~200 K, 
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depending upon various physical parameters such as the infrared emissivity of the comet, 
the penetration of light into the comet's nucleus and the exact composition of the comet 
(Whipple & Huebner 1976). Here, for simplicity, we assume that sublimation from the 
comet follows equation (3) only if its distance from the host star is less than an outer radius, 
Rout', otherwise we assume no sublimation. For numerical purposes, we take R out = 2 AU. 

In writing equation (3), we assume that the sublimation rate from the comet does not 
depend upon whether it is inward-bound or outward-bound. This description is appropriate 
for many Solar System comets (see Whipple & Huebner 1976 and A'Hearn et al. 1995), 
and given the many other unknowns, this seems like the appropriate approximation for the 
models presented here. 

Once the water is ejected from the comet, we assume a radial outflow with a speed, t><> 
For most comets, Vo ~ 1 km s _1 , but for comet Hale-Bopp, Vo was approximately 3 km 
s _1 since the gas was substantially heated as it left the comet's nucleus (Harris et al. 2002). 
Because we are most interested in the unusually massive comets like Hale-Bopp, we use vo 
= 3 km s _1 in our numerical examples. 

If nH 2 oij') denotes the density of water at distance, r, from the comet, then with spher- 
ical symmetry: 

1 d(n H2 or 2 v ) 

~i - - J H 2 on H2 o (4) 

where Jh 2 o(R) denotes the rate at which water is photodissociated. With the boundary 
condition determined by equation (3), the solution to equation (4) is: 

In equation (5), rn 2 o = vh 2 o/ Jh 2 o, is the characteristic distance which the water travels. 
With v = 3 km s" 1 and J/j 2 o(l AU) = 1.0 x 10~ 5 s" 1 (Schleicher & A'Hearn 1988, Harris 
et al. 2002), then r H2 o = 3.0 x 10 10 cm. 

If we define rionij') as the density of OH, then: 
1 d(n 0H r 2 v ) 



r 2 dr 



— —JoH n OH + Jn 2 o n H 2 o (6) 



By writing equation (6), we ignore the relatively small amount of water that is photodisso- 
ciated into O and H 2 . We define rou = vo/Joh, and as long as rou ^ rn 2 o, the solution to 
equation (6) is: 



(7) 
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With v = 3 km s" 1 and J OH (l AU) = 7.5 x 10~ 6 s" 1 (Schleicher & A'Hearn 1988, Harris 
et al. 2002), then tqh = 4 x 10 10 cm. The OH photodissociation rate is a function of solar 
activity and can vary by a factor of 2 (Budzien, Festou & Feldman 1994, van Dishoeck & 
Dalgarno 1984). 

Close to the comet where both r « th 2 o and r « ron, we find from equation (7) 

that: 

Nh 2 o , q x 

uoh ~ -. 8 

A-Kv rH 2 or 

The column density of OH, Noh, along the line of sight with impact parameter b relative to 
the comet's nucleus can be found with the assumptions that uoh is given by equation (8) 
for r < r H . 2 o and it is "small" for r > r H2 o- Thus, for b « r H2 o, 

at N H20 2 r H20 

Noh ~ ^ ln — 7 — 9 

2-kvq ru 2 o b 

Since Nh 2 o varies as R~ 2 and th 2 o varies as R 2 , then Nqh varies approximately as i? -4 . 
Therefore, Noh is sensitive to the comet's orbital distance, R, from the host star. 

We show in Figure 1 the OH column density as a function of displacement from an 
analog to Hale-Bopp for an assumed distance from the Sun of 1 AU and v o = 3 km s -1 . We 
consider both with standard ultraviolet dissociation rates, and du CclSG where Joh and 

Jh 2 o are increased by a factor of 10. In the standard case, the OH cloud is more extended 
and, for small values of b, the column density of OH is relatively low compared to the case 
with high ultraviolet intensity. 

Our results for Noh allow us to determine where the OH is optically thin. As might 
be expected from pumping models for the populatons of the energy levels such as those 
described by Schleicher & A'Hearn (1988), an OH energy level might contain at most 10% of 
all the OH molecules. Then, as seen in figure 1, even in the case of high ultraviolet radiation, 
its maximum column density is less than ~5 x 10 cm" . Since the OH lines near 3100 A 
typically have oscillator strengths near 10~ 3 (Roueff 1996), then with v o = 3 km s" 1 , these 
OH lines are optically thin. However, for R < 1 AU, the OH region is more compact and 
the lines rapidly become optically thick as R decreases. 



2.2. The Comet's Orbital Motion 

The comet's orbit around the star, R(9) ) is given by the expression: 

R(Q) = R p 1 + 6 - (10) 
K ' p 1 + ecos# v ; 
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where R p is the distance at periastron, e is the orbital eccentricity and 6 is the angle measured 
between the major axis of the orbit's ellipse and the host star. The comet's specific angular 
momentum, L, equals RVg where Vq is the tangential speed of the comet. If the mass of the 
star is M*, we can write that: 

L = ([1 + e]GM,R p ) 1 ' 2 (11) 



3. TRANSIENT OH ABSORPTION 

We now compute the absorption produced by the OH cloud around a comet. 



3.1. Threshold of Detectability 

For the conditions of most interest, we expect that the size of the comet's OH cloud is 
smaller than the host star's radius. In this case, we only need to compute the total number 
of OH molecules produced in the outflow from a comet, N(OH). We may use equation (7) 
to find that: 

N(OH) = !n OH A^r 2 dr = ^£ (12) 

J J OH 

Equation (12) shows that N(OH) is independent of the comet's distance to its host star 
because both the water production rate (Nh 2 o) and the OH photodissociation rate {J oh) 
scale as R~ 2 . Equation (12) fails for R > R out where water is not appreciably sublimated. 

The average column density of OH absorbers, N OH in the image of the host star, is 
given by the expression: 

AT - N(OH) 

N OH - (13) 

If OH is optically thin in each portion of the stellar image and if for these exploratory 
models we ignore limb darkening (see Hubeny & Heap 1996), then the total absorption 
in a spectral line just depends upon Noh and not the details on how the gas is spatially 
distributed. Therefore, if the minimum column density for cometary detection is N min , then 
using equations (3), (12) and (13), we require that: 

4n 2 AEJ OH (R)R 2 R 2 N mm 
a > (14) 



Since Jqh varies as R 2 , we see from equation (14) that a is independent of R. 
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The value of N min depends upon the telescopes and instruments that are used. Within 
the interstellar medium where the gas has a similar velocity dispersion to what we expect for 
the outflow from a comet, OH column densities of ~10 13 cm -2 have been measured (Roueff 
1996). Although the OH produced by comets is distributed over more energy levels than 
the OH within the interstellar medium (see, for example, Schleicher & A'Hearn 1988), since 
the same number of photons are absorbed from the background source if the OH is optically 
thin, it may be possible with modern detectors and cross correlation techniques to detect 
total OH column densities in cometary outflows as low as 10 13 cm" 2 . If N min = 10 13 cm -2 , 
we find for environments similar to the solar system's that a > 170 km 2 ; analogs to comet 
Hale-Bopp with a = 1500 km 2 easily satisfy this criterion. Moreover, from equation (14), if 
Nmiu = 10 13 cm -2 , then an analog to comet Hale-Bopp could be detected around a G-type 
star where the ultraviolet luminosity is as much as a factor of 9 larger than the Sun's. 



In our simple model, a comet is only active if R < R ou t- From equation (10), we find 
that for the half of the orbit between 9 = and 6 = ir, the maximum angle for cometary 
activity, 6 max , is: 



For simplicity, we assume that the comet's outflow is undetectable for R < R in where 
Ri n is defined as the region where the OH lines are optically thick. We introduce a minimum 
angle, 8 m i n , in a fashion similar to our definition of 8 max in equation (15) except that we 
substitute i?j n for R ou t- As described below, we are most interested in environments where 
the ultraviolet luminosity of the host star is much larger than the Sun's. Therefore, we take 
R in = 1 AU since, as shown in Figure 1, this is the region where the OH column density 
exceeds 10 14 cm -2 and therefore the outflow from the comet becomes compact and opaque. 
With R in = 1 AU, R out = 2 AU and e xs 1, we list in Table 1 values of 8 min and 9 max for 
analogs to those solar system comets with a > 50 km™ 2 and R p < 2 AU from A'Hearn et al. 
(1995) and also comet Hale-Bopp. 

In order for us to detect transient absorption, the active comet must transit the host 
star when R in < R < R ou t- From the perspective of the host star, we imagine two strips 
of solid angle, one for the inbound and one for the outbound portions of the comet's orbit, 
through which the comet must pass. In each strip, an element of solid angle where the comet 
transits the host star is defined by the angular height, (2R*/R), multiplied by the angular 
width, d6. Thus, the probability that we can witness transient absorption, p, is given by the 



3.2. Likelihood of an Absorption Event 



0, 



max 




(15) 
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fraction of sky covered by the two strips, or, using equation (10): 

2 I" 2 R* R*{[dmax @min} ~\~ e [sin$ ma:E sin$ m j n ]) f-,n\ 

p & — / dd = — — — — (16) 

We list in Table 1 values of p for analogs to recently observed comets; p typically is near 
0.1%. 



3.3. Duration and Probability of Absorption Events 

Above, we have computed the probability that sometime during its orbit a comet can 
be detected. We now estimate the duration, t a b s , of each absorption event. We take for an 
"average" cometary transit time that t a b s = V2R*/Vq. We then write: 

t abs = (17) 
Combining equations (10), (11) and (17), we then find that: 

t m - v / 2^^ /2 ( 1 + ^) 1/2 rif * 

abs[ ] ~ (1 + ecose)(GM*y/i 1 8J 

For e ~ 1, = 7r/2, R p — 1 AU, and M = 1 M star, the characteristic duration of the 
absorption event is 5x 10 4 s. 

The total probability, P(OH), of an observation leading to a detectable transient ab- 
sorption is given by the expression: 

P(OH) = pt abs f (19) 

where T (s^ 1 ) is the rate at which sufficiently large comets arrive with R p < R ou t- Numer- 
ically, with e w 1, 6 = vr/2, R p = 1 AU, M = 1 M and p = lO" 3 , then P(OH) = 47 
t. 

Even for the solar system, T for the large comets that are most likely to be detected 
is uncertain. A'Hearn et al. (1995) identified 4 comets with a > 50 km -2 and R p < 1 
AU during 16 years. For these objects, then T = 8 x 10~ 9 s _1 . In the solar system, the 
cumulative number of comets with area larger than a varies approximately as a -0 7 , at least 
for objects with radii between 1 km and 10 km. (See Meech, Hainaut & Marsden 2004 and 
convert from radius distribution to area distribution). Therefore, scaling from a = 50 km -2 
to a = 1500 km -2 , we expect that comets at least as large as Hale-Bopp have an arrival rate 
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in the inner solar system of 7 x 1CT 10 s~ x . This estimate implies that a very large comet 
like Hale-Bopp should arrive every 50 years, a rate that is approximately consistent with 
the historic data. For example, Sekanina (2002) has argued that comet 1882 Rl had a size 
similar to Hale-Bopp's. With T = 7 x 10~ 10 s _1 , we find from equation (19) that P(OH) = 
3 x 10~ 8 . Thus, the probability of detecting an individual comet around an analog to the 
Sun is very low. 



4. COMETS AROUND YOUNG STARS 

Above, we have considered the possibility of detecting comets around stars similar to 
the Sun. Here, we consider extrapolating these models to models for the Sun at earlier ages. 
There are two competing factors. First, the younger stars are more active and may emit 
as much as 50 times more ultraviolet light (Ayres 1997, Guinan, Ribas & Harper 2003). 
However, the stellar ultraviolet luminosity and level of activity is variable, and some young 
stars have ultraviolet luminosities "only" ~10 times the solar value. Thus analogs to comet 
Hale-Bopp might be detectable. Second, there are both theoretical and observational reasons 
for thinking that young solar-type stars might have vastly greater comet arrival rates than is 
currently witnessed in the solar system. Models for the formation of the Earth all invoke the 
build-up of smaller planetesimals. In this context, Wetherill (1992) has computed that in 
young planetary systems where there is no analog to Jupiter, the rate of cometary periastron 
passages less than 1 AU could be 3 x 10 5 greater than is currently seen in the solar system. 
Also, models for the delivery of water to the Earth describe how giant planets perturb 
the orbits of icy planetesimals into elliptical orbits which result in matter from the outer 
planetary system being injected to within 1 AU of the host star (Morbidelli et al. 2000). 
These planetesimals may be asteroids (Petit, Morbidelli & Chambers 2001) or comets. To 
estimate the magnitude of this effect, consider that the total mass of the Earth's ocean 
is 1.4 x 10 24 g. If this water was mainly delivered during the Earth's first 50 Myr, then 
the water accretion rate onto the Earth was ~10 9 g s _1 . Currently, the zodiacal light is 
explained by a dust infall rate of 3 x 10 6 g s _1 (Fixsen & Dwek 2002) while the Earth only 
accretes about 10 3 g s" 1 of dust (Love & Brownlee 1993). (While dust in the ecliptic plane 
largely results from collisions between asteroids [Grogan, Dermott & Durda 2001], when 
considering the entire sky, an appreciable fraction of the current zodiacal cloud results from 
the disintegration of comets [Hahn et al. 2002]). During the era of ocean formation, if the 
efficiency of the accretion by the Earth of dust with a comet-like origin was 3 x 10 -4 as 
suggested by the current environment, then the inner solar system may have experienced a 
total infall rate of water-rich planetesimals of ~10 12 g s _1 . Therefore, it is possible that in 
the early solar system, the rate of arrival of icy objects within 1 AU was ~3 x 10 5 times 
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larger than the current rate of the arrival of comets. This enhancement of the cometary 
arrival rate suggests that the probability of any single OH-sensitive observation detecting a 
large Hale-Bopp analog in such a young system approaches 1%. 

5. DISCUSSION 

Above, we have discussed models for the detectability of comets. Here, we suggest 
possible targets for an observational program. Although many observers avoid working in 
the near ultraviolet, Israelian, Garcia Lopez & Rebolo (1998) and Boesgaard et al. (1999) 
have observed the OH bands near 3100 A in solar-type stars to infer their oxygen abundances. 
It may be possible to extend such studies to search for comets. 

Comets within 2 AU of a solar-type star not only sublimate ice, but they also eject dust 
into their surroundings. This dust has a characteristic temperature >200 K which results in 
substantial emission at 25 /im. Previous studies with IRAS and ISO of dust around solar-type 
stars have not revealed much infrared excess at A < 25 /im (Laureijs et al. 2002, Zuckerman 
& Song 2004). However, the enormous increase in sensitivity provided by Spitzer (Werner 
et al. 2004) has allowed the study of circumstellar grains around many more solar-type 
stars. Specifically, recent Spitzer observations of 15 Myr old solar-type stars in the Lower 
Centaurus Crux Association show substantial 24 /im excesses where the dust production 
rates may be 10 6 greater than the dust production rate for the zodiacal light within the solar 
system (Chen et al. 2004). These same stars also display relatively low X-ray luminosities 
compared to other solar-type stars in the same Association, and consequently, they may 
also have relatively low ultraviolet luminosities. Therefore they may be suitable targets 
for comet-searching. Some solar- type stars as old as 100 Myr also show substantial 24 /im 
excesses (Gorlova et al. 2004). 

Unfortunately, the prospects for detecting transient OH absorption around warmer stars 
like (3 Pic (L = 8.5 L ; T eff = 8200 K; Crifo et al. 1997) are not very good. That is, (3 
Pic emits enough photospheric light between 1600 A and 1800 A, where both OH and water 
have significant photodissociation cross sections, that the lifetimes of these molecules are 
relatively short. In particular, given the star's luminosity, we would expect that comets at 3 
AU around (3 Pic might resemble comets at 1 AU around the Sun. Using archived HST data 
(see Roberge et al. 2000), scaling the observed ultraviolet flux upwards by (19.3 pc/3 AU) 2 , 
and using the compilation of the photodestruction cross sections by Budzien et al. (1994), 
we find that the lifetimes at 3 AU around (3 Pic of OH and H2O are ~1600 s and ~330 s, 
respectively. As a result, from equation (14), only comets with areas greater than 3 x 10 4 
km 2 would be large enough to produce an OH column density of 10 13 cm -2 . Thus transient 
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OH absorption spectrum is expected to be very weak unless the comet is vastly larger than 
Hale-Bopp. Some gas-phase CO is seen in the ultraviolet absorption- line spectrum of (3 Pic, 
but it does not change with time and it is probably not directly associated with infalling 
comets (Roberge et al. 2000). 

6. CONCLUSIONS 

We have considered the transient OH absorption signature of comets similar to those 
found in the solar system. Such events are very rare for stars like the Sun. For young solar- 
type stars with a substantial 24 fim excess where there may be comets delivering water to 
terrestrial planets, a single, sensitive, randomly-timed observation may have a 1% chance of 
detecting OH flowing away from a large comet. 

This work has been partly supported by NASA. 
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Table 1 - OH Detectability of Analogs to Recent Bright Comets 



Comet 


R p 


a 




@max 


P 




(AU) 


(km 2 ) 


(°) 


(°) 




Bradfield 


0.26 


96 


119 


138 


0.00036 


Halley 


0.59 


50 


80 


114 


0.00066 


Levy 


0.94 


110 


28 


93 


0.0013 


Meier 


1.14 


250 





82 


0.0016 


Wilson 


1.20 


79 





78 


0.0015 


West 


0.20 


130 


127 


143 


0.00031 


Hale-Bopp 


0.91 


1500 


35 


95 


0.0012 



We take max and p from equations (15) and (16), respectively. The cometary parameters 
are from A'Hearn et al. (1995) except for Hale-Bopp which is discussed in the text. 
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Fig. 1. — Plot of OH column density vs. impact parameter for an analog to comet Hale-Bopp 
with a = 1500 km 2 located at 1 AU from a star with L — 1 L . We use vo — 3 km s^ 1 both 
for the case with standard photodissociation rates (solid line) and for the case where these 
rates are increased by a factor of 10 (dashed line). 



